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Abstract 

Motivated by a recent experiment that reported the successful synthesis of hexagonal (/i) AIN 
[Tsipas et al. Appl. Phys. Lett. 103 , 251605 (2013)] we investigate structural, electronic and 
vibrational properties of bulk, bilayer and monolayer structures of /i-AlN by using first-principles 
calculations. We show that the hexagonal phase of the bulk h-AlN is a stable direct-bandgap 
semiconductor. Calculated phonon spectrum displays a rigid-layer shear mode at 274 cm“^ and 
an Eg mode at 703 cm“^ which are observable by Raman measurements. In addition, single layer 
/i-AlN is an indirect-bandgap semiconductor with a nonmagnetic ground state. For the bilayer 
structure, A A' type stacking is found to be the most favorable one and interlayer interaction is 
strong. While N-layered /i-AlN is an indirect bandgap semiconductor for N=l-10, we predict that 
thicker structures (N>10) have a direct-bandgap at the F-point. The number-of-layer-dependent 
bandgap transitions in h-AlN is interesting in that it is significantly different from the indirect-to- 
direct crossover obtained in the transition metal dichalcogenides. 

PACS numbers: 81.05.ue, 85.12.de, 68.47.Fg, 68.43.Bc, 68.43.Fg 
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I. INTRODUCTION 


Bulk structures of III-V semiconductors have been widely studied due to their impor¬ 
tance for technological applications such as blue light-emitting diodes, lasers operating in 
the blue and the ultraviolet regime, high temperature diodes and transistors.®^ In addi¬ 
tion, technological advances lead to the emergence of novel low-dimensional forms of III- 
V binary compounds. Experimental fabrication of AIN nanowires,®^ nanobelt^BII^ and 
nanodot^^^nSl have already been reported. In recent years, theoretical and experimental 
studies of graphend^^ provided a wide range of knowledge for a new class of materials and 
have opened up possibilities for the synthesis of many similar structures such as; silicene,^^®^ 
germanene,®® transition metal dichalcogenides (TMDs)P^^ and hexagonal structures of 
III-V binary compounds (e.g. h-BN, h-AlN).®®^ Among the binary compounds, after single 
layer BN, AIN is currently a material of much interest due to its semiconducting nature and 
its large bandgap which is suitable for device applications. 

Following the hrst theoretical study on the stability of hexagonal AIN reported by Sahin 
et successful experimental realization of hexagonal AIN phase was achieved very re¬ 
cently by Tsipas et alW^ In the early study of Du et ai, energetic and electronic properties 
of one-dimensional AIN nanostructures such as nanowires with hexagonal cross sections, 
double-and triple-walled faceted nanotubes and single-walled faceted AIN nanotubes were 
investigated.!^ Between the hrst theoretical prediction and the experimental realization of 
h-AlN, other groups have also focused on this material. Zheng et a/.®! predicted that zigzag 
AIN nanoribbons have an indirect bandgap whereas armchair AIN nanoribbons have a di¬ 
rect bandgap, and these bandgaps monotonically decrease with increasing ribbon width. 
Almeida et al^ studied the energetics and electronic properties of typical defects in a h- 
AIN network such as vacancies, anti-sites and impurities. It was shown that defects such as 
N vacancies and Si impurities lead to the breaking of the planar symmetry of the h-AlN sheet 
and signihcant changes in the band structure in the vicinity of the Fermi level. In addition, 
Chen et al. systematically investigated the electronic structure of armchair and zigzag AIN, 
GaN nanoribbons and investigated also the electronic properties of AlN/GaN nanoribbon 
heterojunctions. They found that the bandgap of both nanoribbons decrease monotonically 
as the ribbon widths increase and that the bandgap of the nanoribbon heterojunctions are 
closely related to the AlN/GaN ratio.^^ Shi et al. calculated the magnetic properties of 
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FIG. 1: (Color online) Top and tilted-side views of atomic structures of bulk (a) wurtzite and (b) 
hexagonal AIN (top panel). Their electronic band dispersion (lower panel). 

undoped and transition metal (TM) doped AIN nanosheets by using first-principles calcu¬ 
lations. They reported that AIN nanosheet is nonmagnetic, whereas a single 3d TM atom 
can bring about large local magnetic moments in TM-doped AIN nanosheets.^ 

Motivated by the recent study of Tsipas et alW^ reporting the formation of stable h- 
AIN phases at the early stages of AIN growth, we aim to address the following experimental 
observations: (i) modihcation in electronic valence band structure, (ii) relatively large lattice 
parameter that was reported as an indication of h phase, (hi) reduction in band gap in 
comparison with bulk wurtzite, and (iv) dynamical stability of h-AlN. 

Our study reveal that the synthesized h-AlN, although it has structural similarities, 
electronically displays different characteristics from h-BN and TMDs. As well as weak 
van der Waals (vdW) interactions, ionic character of the interlayer interactions plays an 
important role in the electronic properties of mutilayered h-AlN. Therefore, while monolayers 
of h-BN and TMDs are direct band gap semiconductors, h-AlN has an indirect band gap. 
In addition, upon dimensional reduction from bulk to bilayer energy the band gap decreases 
like as was found for h-BN and TMDs. 
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The paper is organized as follows: in Sec. [TTjwe give details of onr compntational method¬ 


ology. In Sec. Ill we briefly overview the strnctnral and electronic properties of hexagonal 


AIN (h-AlN) strnctnre together with wnrtzite strnctnre. In Sec. IV stability of few layer 
strnctnres and their electronic properties are investigated in details. Onr findings are dis- 
cnssed in Sec. El 


II. COMPUTATIONAL METHODOLOGY 

Onr investigations of the minimnm energy confignrations of different AIN strnctnres and 
their electronic properties, were carried ont nsing the Vienna ab-initio simnlation package 
VASP.H3H1ZI which is based on density fnnctional theory (DFT). The VAST code solves the 
Kohn-Sham eqnationP^ for a system with periodic bonndary conditions nsing plane-wave 
basis set iteratively. The Perdew-Bnrke-Ernzerhof (PBE) form of the Generalized Gradient 
Approximation (GGA)3^^ was adopted to describe electron exchange and correlation. The 
band gap nnderestimation in bare GGA calcnlations are also examined by nsing the hybrid 
DFT-IISE06 fnnctional.^ The vdW interaction that is significant for layered h-AlN material 
was taken into acconnt for all mnltilayer strnctnres.l^^^ For the charge transfer analysis the 
Bader techniqne was nsed.^^ 

The kinetic energy cnt-off of the plane-wave basis set was 500 eV in all calcnlations. 
The optimization of atomic positions was performed by minimizing the total energy and 
the forces on the atoms hence the energy difference between seqnential steps was taken 
10“^ eV as a convergence criterion in the strnctnral relaxation, and the convergence for the 
Hellmann-Feynman forces on each atom was taken to be 0.05 eV/A. In addition, Ganssian 
smearing factor of 0.05 eV was nsed for non-self-consistent calculations and the pressures 
on the unit cell were decreased to a value less then 1.0 kB in all three directions. For the 
determination of accurate charge densities, Brillouin zone integration was performed using a 
35 X 35 X 1 Monkhorst-Pack mesh for the primitive unit cell.^ To avoid interaction between 
adjacent AIN monolayers and few layer systems, our calculations were performed with a 
large unit cell including 14 A vacuum space. We also calculated the cohesive energy 
(Ecoh) which was formulated as E^oh = {nAiEAi + unEj^ - E^) / {hai + un), where E^, Eai, 
Eat, hai, and are the total energy per unitcell, the energy of free A1 atom, the energy of 
free N atom, number of A1 and N atoms in unitcell, respectively. 


4 








III. STRUCTURES OF BULK ALN 


First we briefly overview the characteristics of wurtzite and recently synthesized h-AlN 
before the extensive investigation of the number of layer dependent properties of h-AlN 
are undertaken. At ambient conditions, AIN crystallizes in a wurtzite {wz) structure that 
belongs to the PQ^mc space group. As shown in Table [T| PBE approximation gives the 
following structural parameters a = 3.11 A and ^ = 1.61. In the wz structure of AIN 
each Al-N bond is formed by 2.35 e~ charge donation from A1 to N atom and therefore 
it has a highly ionic character. In parallel with the abundance of the wz form of AIN in 
nature, among the possible bulk structures wz is the energetically most favorable one. The 
electronic dispersion shown in Fig. indicates that tcz-AlN is a semiconductor with a 4.2 
eV direct bandgap at the F point. 

TABLE I: Calculated lattice parameter in the lateral direction a and lattice parameter in the 
vertical direction c, the distance between layers d^i, the intralayer atomic distance between A1 
and N atoms dAi-N, the charge transfer from A1 atom to N atom Ap, the cohesive energy E^oh 
and the energy bandgap of the structure Egap- (d) and (i) indicate direct and indirect bandgap, 
respectively. 


a c d-LL d^i-N 

(A) (A) (A) (A) 

^coh ^gap 

(e) (eV) (eV) 

Bulk tcz-AlN 3.11 5.01 

1.90 

2.35 6.14 4.2 (d) 

Bulk h-AlN 3.30 4.15 2.08 

1.90 

2.37 6.02 3.4 (d) 

2L h-AlN 3.20 - 2.13 

1.85 

2.32 5.73 3.5 (i) 

IL h-AlN 3.13 - 

1.81 

2.28 5.36 2.9 (i) 


In addition to some early theoretical predictions^ the stability of the hexagonal phase of 
AIN (h-AlN) was proven by a very recent experimental study.^^ Regarding the less layered 
structures or surfaces of wurtzite materials, transformation from wurtzite to a graphite-like 
structure that allows the removal of destabilizing dipoles is energetically more favorable.^ 
This graphite-like hexagonal structure of AIN belongs to the Pd^/mmc space group. The lay¬ 
ered planar structure has two atoms in the unitcell which has lattice vectors vi = a(^, 0 ), 

V 2 = a(|, —^, 0 ) and V 3 = c( 0 , 0 , 1 ) where |vi| = |v 2 |. Atomic coordinates are (^, ^, 0 ) 
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and ^^,0), for first and second type of atoms respectively. We considered that the 

layers are AA' stacking (deciding the stacking of layers is explained in Section IV-B) which 
is shown in Figj^b). And we found that the lattice constant is a = 3.30 A and the distance 
between layers is c = 2.08 A which is smaller when compared to graphite and h-BN (~ 3.33 
A and ~ 3.30 A for graphitd^*^ and respectively). Intraleyer atomic distance 

between A1 and N atoms is 1.90 A. This is the same as for tc^-AlN. The cohesive energies 
are given in Table 

As shown in Fig. [^d), layered h-AlN is a direct bandgap semiconductor which has 3.4 
eV bandgap with band edges located at the F point. Therefore, our results conhrm the 
experimentally reported bandgap reductiorP^ in h-AlN compared to bulk taz-AlN. At the 
top of the VB, there is a single parabolic band for holes and the main contributions are from 
the Pz orbital of N atoms. The bottom of the CB, there is a single parabolic band and most 
of the states come from p^. orbital of N atoms. Final charges of A1 and N atoms are 0.63 
e~ and 7.37 e~, respectively. Therefore, 2.37 e~ charges are transferred from A1 to N atom. 
Following section is devoted to the understanding of the characteristics of monolayer h-AlN 
and the nature of the inter-layer interactions. 



FIG. 2: (Color online) Phonon dispersion of bulk h-AlN. 

Lastly, we analyze vibrational properties of the h-AlN phase by calculating the phonon 
spectrum in whole Brillouin Zone (BZ) using the small displacement methodologjP^ for 
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FIG. 3: (Color online) Monolayer /i-AlN band diagram within PBE and HSE06 (left panel), 
partial density of states (middle panel) and the band decomposed charge densities of the band 
edges (right panel). Sections of 3D charge density that represent the connection to the neighboring 
cells are filled with blue color. 

32-atomic supercell. As shown in Fig. [^phonon spectrum of h-AlN displays real eigenfre- 
quencies in the whole BZ. Therefore, there is no doubt on the dynamical stability of the 
h-AlN structure. Here we also see that differing from its h-BN counterpart the optical and 
acoustical modes couple with each other at certain points in the BZ. In the case of h-AlN, 
there are 4 atoms per primitive unit-cell and 12 phonon branches. The hrst three are acous¬ 
tical phonon modes, the phonon branch at 274 cm“^ is a doubly degenerate rigid layer shear 
mode, known as low-energy mode, involving the out-of-phase motion of atoms in adjacent 
planes. In addition, the highly dispersive phonon modes with eigenvalues 406 and 679 cm“^ 
at the zone center correspond to out-of-phase motion of A1 and N atoms in adjacent layers. 
Interestingly, the contribution of the second atom type in these modes are negligibly small. 
The highest optical mode, doubly degenerate at the T point, have the symmetry and 
therefore it is expected to be measured in Raman experiments. 
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IV. MONOLAYER, BILAYER AND FEW LAYER /i-AlN 


A. Monolayer /i-AlN 

The first prediction of the dynamical stability and electronic properties of single layers of 
h-AlN and similar III-V compounds were hrst reported by Sahin et a/® Monolayer hexagonal 
structure belongs to the space group PQ^/mmc with unit vectors vi = a(|,^,0), V 2 = 
a(|,—^,0) where |vi| = |v 2 |. In this conhguration the atomic coordinates are given as 
0) and 0) for first and second type of atoms, respectively. We calculated 

that the lattice constant of monolayer h-AlN is 3.13 A and the distance between A1 and 
N atoms is 1.81 A which is the lowest value when compared with bulk forms and bilayer 
structure as seen from Table [1} These results are in good agreement with the measured lattice 
parameters (3.13 A) of epitaxially grown of h-AlN, at early stages, by Tsipas et a/..^The 
cohesive energy of monolayer h-AlN is the lowest among the possible phases. It is also seen 
that the amount of charge transfer from A1 to N slightly decreases from bulk to monolayer 
h-AlN. 

Fig. i illustrates the band diagrams, density of states and also the 3D charge densities 
of states at the band edges for monolayer h-AlN. Differing from bulk h-AlN which has a 
direct bandgap at the T, the VBM of the monolayer is at the K point and the CBM is at the 
T point, therefore monolayer h-AlN is an indirect bandgap insulator with 2.9 eV bandgap. 
It is also seen that bare-GGA and DFT-HSE06 approximated electronic band dispersions 
shown in Fig. are almost the same. With DFT-HSE06 hybrid functional, h-AlN is an 
indirect bandgap material with a bandgap of 4.06 eV. It is also noteworthy to mention that 
h-BN is the counterpart of AIN with a direct bandgap at the K point. 

Right panel of Figshows that the bonding states at VBM (label I) are mainly composed 
of N-p^ orbitals and the degenerate state, with lower energy, at the F point (label II) stems 
from hybridized Px and Py orbitals. In addition, the CBM located at the F point has a 
quadratic-like dispersion (label III) and the relevant charge density is in the form of surface 
state. Therefore, the conduction electrons in single layer h-AlN display free-electron-like 
behavior. Moreover the anti-bonding states with higher energy at the M and the K points 
are composed by mainly A1 and slightly N states together. 
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FIG. 4: (Color online) (a) Possible stackings between two layers of /i-AlN. (b) Layer-layer in¬ 
teraction energy for different stackings, (c) Interlayer spacing dependent band dispersion of A A' 
stacked two layers of /i-AlN. 

B. Bilayer /i-AlN and Layer-layer Interaction 

During the growth of lamellar materials prediction of the most preferable stacking se¬ 
quence is of importance. To determine the stacking order of bilayer AIN, we examined hve 
possible (see Fig. |^a)) stacking types: AA (A1 over A1 and N over N), AA' (A1 over N and 
N over Al), AB (A1 over N and N over center of hexagon), A'B (A1 over A1 and N over center 
of hexagon), AB' (N over N and Al over center of hexagon). 

To determine the interaction strength and interlayer interaction prohle of above- 
mentioned stacking orders of bilayer h-AlN, total energies were calculated as a function 
of layer-layer distance. As seen from Fig. |^b), the largest interlayer coupling energy and 
shortest interlayer distance is obtained for AA' type stacking. Therefore AA' stacking cor- 
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FIG. 5; (Color online) Evolution of electronic band dispersion of N-layer /i-AlN. 

responds to the ground state stacking order. It is also seen that the vdW correction to 
the interaction energy is ~250 meV per unitcell. In addition, for AA' stacking, layer-layer 
distance is 2.13 A which is the minimum value among all stacking combinations. Here it 
appears that due to the presence of surface states interlayer distance in bilayer h-AlN is 
slightly larger than that of bulk h-AlN. However, compared to graphit^®^ and 
and their bilayer forms, h-AlN has a smaller interlayer distance. The lattice parameter, 
3.20 A, which is the maximum value of all stacking combinations that slightly differs from 
the value obtained by self-correlation analysis of the experimental data (3.14±0.06) for few 
layers. This slight disagreement can be understood from the effect of the Ag(lll) substrate. 

Moreover, in Fig. |^c), we analyze how the electronic band dispersion of AA'-stacked 
bilayer is modihed with varying interlayer distance. Negligibly interacting, weakly interact¬ 
ing and highest interaction (ground state) cases are illustrated using layer-layer distances of 
5.00, 3.00 and 2.13 A, respectively. The band diagram of negligibly interacting case resem¬ 
bles the electronic structure of monolayer h-AlN. However, when the layers start to interact, 
for instance in the weakly interacting case (delineated by green lines in Fig. |^c)), the top 
of the VBM are shifted towards the F zone center. When the interlayer distance reaches its 
optimum position, one of the VBM appears at K-F while the other one is located at M-F. 
Interestingly, the location of CBM in BZ is independent from the interlayer distance. Hence 
from monolayer to bilayer, the indirect character is conserved and the bandgap is changed 
from 2.9 to 3.5 eV. 
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TABLE II: Calculated direct and indirect transition energies between the VB and CB edges of 
single to ten layered structures. 



r —> 

r K—^ 

r TM—^ 

r KF —^ 

F Eco/i 

IL AIN 

3.62 

2.92 

- 

- 

5.36 

2L AIN 

4.29 

- 

3.61 

3.53 

5.73 

3L AIN 

4.05 

3.58 

3.57 

3.53 

5.82 

4L AIN 

3.79 

- 

3.51 

3.50 

5.86 

5L AIN 

3.61 

3.73 

3.46 

3.46 

5.89 

6 L AIN 

3.48 

- 

3.41 

3.41 

5.92 

7L AIN 

3.38 

3.77 

3.35 

3.35 

5.93 

8 L AIN 

3.33 

- 

3.30 

3.30 

5.94 

9L AIN 

3.27 

3.81 

3.26 

3.26 

5.95 

lOL AIN 

3.22 

- 

- 

- 

5.96 


C. Few Layer /i-AlN 

Synthesis of stable ultra-thin h-AlN structures, sub-monolayer to 12 layers, was demon¬ 
strated by Tsipas et a/. recently. Our first-principles total energy optimization calculations 
revealed that lamellar hexagonal structure of N-layered (N=l,2,...) AIN is a stable phase. 
As shown in Figs. [^and|^ planarity of the few-layered structures is preserved, except the 
negligible buckling at the surface layers, while electronically few layer h-AlN is significantly 
different from monolayer one. 

In order to illustrate how the electronic structure evolves with increasing number of layers 
we present the band dispersion of 3, 5, 8 and 10-layered h-AlN in Fig. Here there are 
several interesting trends that is worth to note: (i) the energy of the conduction band edge 
at the K point decreases with increasing number of layers, (ii) states at conduction band 
edge which have surface state character remain unchanged. The energy difference between 
the band edges VBM(iF) and CBM(F) increases from 2.92 eV for monolayer to 3.81 eV 
for 10-layer h-AlN, and (iii) the most significant influence of additional number of layers is 
on the bonding states forming the valence band maximum at the FM and the FK points. 
It appears that due to the strong interlayer interaction each additional layer shifts these 
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FIG. 6: (Color online) The evolution of the bonding charge density of /i-AlN for monolayer, 
bilayer, 3-layer, 8-layer and 15-layer. Sections of 3D charge density that represent the connection 
to the neighboring cells are filled with blue color. 

band edges towards the zone center. This gradnal shift of VB edges continnes np to 10 
layers. Upon the formation of the 10th layer h-AlN strnctnre reaches the bnlk limit and 
it becomes a direct band gap semicondnctor. Increasing the thickness fnrthermore has no 
influence on the electronic properties of N-layered h-AlN and they are direct bandgap (3.22 
eV for ten layered) materials like bnlk h-AlN. As given in detail in Table transition 
energies between the points converge asymptotically to valnes of bnlk h-AlN. Increasing the 
nnmber of layers, cohesive energies converge to that of bnlk h-AlN. It appears that only 
few layer h-AlN strnctnres with thickness N>11 are snitable for lasing device applications. 
Here calcnlated modifications in electronic strnctnre explains and snpports the experimental 
findings of Tsipas et al.W^ 

However, the evolntion of the VBM with changing nnmber of layers reqnires farther at¬ 
tention. In Fig. [^we present the 3D charge density of the electronic states located at the 
VBM. It is clearly seen that the character of the bonding state (VBM) is modified monoton- 
ically with increasing nnmber of layers. With increasing nnmber of layers the hybridization 
between N-p^ states is weakened. In addition, dne to the mild bnckling the interaction of 
the snrface states which have sp^-like character and states vanishes at the vicinity of 
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the upper- and the lower-most layers. Therefore for enough thick h-AlN materials, which is 
found to be larger than 10-layers, the electronic structure is mainly determined by uniformly 
distributed N-p^ states and hence the structure behaves like a bulk material having a direct 
bandgap at the T symmetry point. 

V. CONCLUSIONS 

Using hrst principles methodology, we investigated the thickness dependent electronic 
properties of layered hexagonal AIN. First, we showed that the bulk structure of hexagonal 
AIN is a semiconductor with a direct bandgap at the F point. Phonon spectrum analysis 
of this structure reveals two Raman-active modes at 274 and 703 cm“^. Here, the lattice 
parameter of h-AlN, which is larger than the wurtzite phase, agrees with the experiment 
of Tsipas et that reports the formation of hexagonal regions at the early stages of the 
growth process. Next, the formation of multilayer structures and most favorable stacking 
order were investigated via total energy calculations. It is seen that, similar to the hexagonal 
BN counterpart, among the possible stackings AA' type stacking is the most favorable. It is 
worthwhile to note that the interlayer interaction is larger compared to similar layered ma¬ 
terials with hexagonal lattice symmetry. Differing from some other two-dimensional crystal 
structures such as graphite, h-BN and TMDs, here the ionic interlayer interaction is domi¬ 
nant between the h-AlN layers. Therefore, for the synthesis of h-AlN the epitaxial growth 
technique performed by Tsipas et appears more suitable than mechanical exfoliation. 

Subsequently, the evolution of the electronic structure of N-layered h-AlN was investi¬ 
gated for structures with 1-15 layers. It is seen that unlike bulk h-AlN, monolayer h-AlN 
has an indirect bandgap where VBM and CBM are at the K and the F points, respec¬ 
tively. Moreover, it is seen that upon the formation of additional layers valence band edges 
gradually shift towards the Brillouin Zone center. Such a modihcation of the valence band 
states is in good agreement with the experiment of Tsipas et It was also calculated 
that the few layered structures reach the bulk limit and their electronic properties remains 
almost unchanged after the formation of 10 layers. Therefore, structures thicker than 10 
layers, exhibits a direct band gap of 3.22 eV at the F point. We believe that easy synthesis 
and the presence of a thickness-dependent bandgap crossover in few-layered hexagonal AIN 
structures are very important for novel device applications. 
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